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PART  ONi:.  EXACT  DETECTION  I’ROUAlilLlTY 


I.  INTRODUCTION 


The  determination  of  detection  probability  when  the  sum  of  N  detected 

pulses  of  signal  plus  noise  is  compared  to  a  threshold  hi.s  been  studied  by 

Marcum,  Swerling,  Schwartz,  and  Vanntcola,  among  others.  Marcum^  treated 

2 

the  non- fluctuating  target;  Swerling  treated  two  limiting  fluctuating  tar¬ 
get  situations  (complete  correlation  and  complete  decorrelation)  -  he  also 
presented  a  general  method  for  treating  arbitrary  correlation^,  but  aave  no 
results  for  particular  correlation  models;  Schwartz  achieved  an  exact  solu¬ 
tion  when  there  are  only  two  pulses;  Vannicola^  has  constructed  a  solution 
from  Swerling' s  two  extreme  fluctuation  models  by  considering  M  independent 
sets  of  N  fully  correlated  pulses  -  his  solution  applies  only  to  targets 
which  are  block  correlated.  This  paper  extends  the  previously  referenced 
work  and  presents  exact  results  for  a  Rayleigh  target  whose  inphase  (and 
quadrature)  components  have  exponential  correlation. 


1 >2,3,4, References  at  end  of 


text . 


1 


11.  ANALYSIS 

bcrote  the  inphase  and  quadrature  components  of  signal  and  noise  respect¬ 
ively  by  the  Nxl  complex  vectors  A+iB  and  X+iY.  We  assume  (i)  A,D  are  i.i.d.; 
(ii)  X,Y  are  i.i.d  stationary  Gaussian;  (iii)  the  noise  is  additive  and  inde¬ 
pendent  of  the  signal.  Employing  the  Noyman-Pearson  detection  criterion,  the 
probability  of  detection  is  given  by 


A  ^ 

P_  -  J  p  (v)dv 

D  •'  V 


where  v  is  the  integrator  output  normalized  to  the  average  per  pulse  noise 
power ,  i . e . , 


A  Ia+xI^  +  Ib+yI^ 


and  V  is  the  normalized  thresho' i  determined  by  the  specified  probability 
T 

of  false  alarm. 


10  I 

P  ^  f  9  (v Ma!^  +  Ib'i^  -  0)  dv 


In  order  to  determine  we  require  the  probability  density  function  (pdf) 
of  V.  We  first  write  P^(v)  as  an  inverse  Laplc -e  transform 


sv 

L  (s)e  ds 

V 


U 


2 


Next  recognizing 


■  / 


p^(v)e  *'^dv 


(5) 


~*8V 

as  the  expectation  of  e  we  proceed  to  calculate  this  expectation  over 
the  domain  of  A,B>X,Y  i.e.. 


A+X  +  B+Y 


-////n(A,B,X,Y)e-®  ^^2 


dAdBdXdY 


(6) 


In  view  of  the  three  assumptions  characterizing  the  signal  and  noise,  the 
joint  pdf  of  A,B,X,Y  may  be  written  as 

Ix]^  +  IyI^ 


d(A,B,X,V)  -  p. (A)p. (B) 

■  rt  n 
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<2^6"") 


(7) 


Thus  (6)  becomes 


L  (s) 

V 


«.00  ^OJ 


(s1a+x1^  +  !v|2) 

26 


(2tt6^) 


n72 


dX 


(8) 


or,  completing  the  square  in  X, 


L^(s)  - 


00 

/ 


s+1 


Pj^(A)e 


n 

26‘ 


-O)  (s  +  l) 


N/2 


dA 


/ 


S+l  I „  s  ,  I  2 
[2tt6^/(s+1)]’^''“ 


dX 


(9) 


I  *  . 


w  •  - 


^  1 


Thus  for  any  signal  model  we  have  the  general  formula 


L^(s) 


CD 

-.1. 


(s+1)  L-s» 


p^(A)e  dA 


A.  Marcum  Model 


The  non-fluctuating  target  may  be  accommodated  by  choosing 


. V  “ 


where  6  (•)  is  the  Dirac  delta  function.  Then  (10)  yields  Marcum's 
result 


,  ,  ,  1  s+1 

(s+1) 


where  x  is  the  per  pulse  ratio  of  signal  power  to  average  noise  power,  i.e.. 


^  “  “2 
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B.  Sworl ing  Models 

Typical  radar  targets  are  composed  of  a  large  number  of  scatterers 
whose  distances  relative  to  the  radar  change  with  body  vibration  and  vary¬ 
ing  target  aspect.  The  echoes  from  the  individual  scatters  contribute 


.*  .•w  ,*«  *.  .*•  V  N.  ♦. 


,4 


f 


t 


« 


I 


amplitude  and  phase  terms  which  combine  at  the  radar  frequency  to  produce 
a  fluctuating  signal.  Swerling  has  bounded  the  effects  of  target  correla¬ 
tion  by  considering  two  limiting  situations:  The  alow  fluctuation  model 
(unity  air4)litude  correlation  on  a  scan,  scan-to-scan  independence)  and  the 
fast  fluctuation  model  (zero  amplitude  correlation  on  a  scan,  pulae-to- 
pulse  independence) .  These  may  be  accommodated  by  choosing  respectively 

N 

p^  (a^,...,a^)  “  (14) 

and 


Pa  <“l . V  “  I  Pa  V  '15) 

6 

Measurements  of  echoes  from  ;aany  types  of  radar  targets  confirm  that  the 
best  characterization  of  the  first  order  amplitude  statistics  is  given  Ly 
the  Rayleigh  distribution.  (The  best  characterization  of  the  autocorrela¬ 
tion  function  of  the  radar  cross  section  is  given  by  an  exponential  func¬ 
tion.)  Thus  we  assume 


P  _ 

Va^b 


(r) 


(0  <  r  «>) 


(16) 


where  20*'  is  the  average  target  cross  section.  In  the  remainder  of  this 
paper  we  confine  our  attention  to  this  "Rayleigh  target". 

The  inphase  (or  quadrature)  components  have  first  order  statistics 
given  by 


P^(a, 


(17) 


See  Roforeucos  at  end  of  text. 
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Thus  the  general  formula  (IQ)  specializes  to  Swerling's  slow  fluctua¬ 
tion  result  (case  1) 


L  (3) 
V 


(s+1)**  ^  ((l+NX)s+l] 


when  (14)  is  employed  as  the  correlation  model  and  to  Swerling's  fast 
fluctuation  result  (case  2) 


L^(s) 


t(l+x)s+l]‘ 


when  (15)  is  employed  as  the  correlation  model. 


C.  Partially  Correlated  Model 

In  order  to  accommodate  Rayleigh  targets  which  give  rite  to  partially 
correlated  signals  (0<p<l)  we  introduce  the  general  correlation  model 


2a 

. *»’  ■  lAi* 


where  "T"  denotes  transpose  and  C  is  the  covariance  matrix  of  t 
For  this  correlated  Rayleigh  model,  (10)  yields 


Swerling's  results  may  be  obtained  from  this  formula  by  choosing  ^^^”0 

(nj*!)  for  the  slow  fluctuation  case  [c.f.(18)],  and  X  =1  for  the  fast  fluctu- 

n 

ation  case  [c.f-ClS)]-  These  correspond  respectively  to  the  singular  co- 
variance  matrix 


C 


P-1 


(25) 


in  which  the  first  pulse  contains  all  the  averai^e  signal  power  of  the 
pulse  train  and  the  remaining  pulses  contain  no  signal,  and  to  the  identity 
matrix 


C 


p«0 


(26) 


in  which  the  total  average  signal  power  is  liniformly  distributed  eunong  the 
N  pulses. 

The  inverse  r.aplace  transfor.m  of  (24)  now  yields  the  pdf  of  the  inte¬ 
grator  output.  When  X”0  we  obtain 


p^(vlx»0) 


-V 

e 


N-l 

V _ 

(N-l) ! 


and  when  (assuming  distinct  eigenvalues)  we  obtain 


(27) 


8 


V 


(28) 


Thus  (c,f.(3)]  the  threshold,  V  ,  is  given  in  terms  of  the  specified  P 

T  FA 

by 


P 


FA 


ll“U 


n  t 


(29) 


and  [c.f.(l)l  the  probability  of  detection  by 


-D-r 


U«1 


e 


l^xXn 


(30) 


In  order  to  calcula  e  for  any  particular  correlation  model  we  must 
provide  a  description  of  the  correlation  matrix.  If  the  signal  arises  from 
a  stationary  process,  then  C  is  a  symmetric  Toeplitz  matrix  with  N  distinct 
elements.  We  should  IDte  to  associate  the  eigenvalues  witJi  a  single  correl¬ 
ation  parameter,  p,  in  terms  of  which  may  be  conveniently  characterized. 

For  this  purpose  we  assume  that  the  signal  is  described  by  a  first  order 
Markov  process.  This  is  consistent  with  Edrington's  measurements  which  show 
that  the  radar  cross  section  is  exponentially  correlated. 

Consider  a  signal  which  consists  of  a  train  of  N  pulses  with  uniform  spac¬ 
ing  T.  Let  the  k,n  element  of  C  be  taken  as 

_  -  1 k-n I  Tv  A  j  k-n | 

C  *  e  ^0 

kn 


(0<P<1) 


(31) 


than 


C  - 


f  ^  I 


N-l 


(32) 


The  eigenvalues  of  C  provide  a  non-trivial  solution  to  the  matrix  equation 

tC  -  Xllu  -  0  (33) 

Since  the  sum  of  the  eigenvalues  equals  the  trace  of  C  we  have  the  relation 


N 

y  X  -  N 

,  n 


(34) 


which  we  make  use  of  later. 

Consider  first  the  special  case  of  a  low  prf  (pulse  repetition  frequency) 
waveform  for  which  the  interpulse  spacing,  T,  is  so  large  that  the  correla¬ 
tion  of  noncontiguous  pulses  may  be  neglected.  Thus  C  becomes  tri-diagonal 
and  (33)  is  equivalent  to  a  homogeneous  bvp  (boundary  value  problem) .  This 
consists  of  a  set  of  homogeneous  second  order  difference  equations 

pu  ,  +  (l-X)u  +  pu  ,  ■  0  (n»l,...,N)  (35) 

n-l  n  n+1 


together  with  the  homogeneous  boundary  conditions 


0 


(36) 


Since  the  equation  is  linear  and  has  constant  coefficients,  there  are 


10 


two  solutions  in  the  form 


where 


<37) 


y 


X-1 

2P 


± 


(38) 


The  condition 


A-i 


2p 


^  1  implies  either  X  ^  l+2p  or  X  <  l-2p,  each  of  which 

X-ll 


leads  to  a  contradiction  of  (34) .  Hence  we  must  have 
±i6 

Y  ■  e  where 


2p 


<  1,  i.e. , 


A  ^ 
”  2P 


(39) 


We  write  the  general  solution  to  (35)  in  terms  of  e  as 

u  “  K  cos  n6  +  K  sin  n9  (40) 


Since  Uq  *  0,  we  have 


u  “  K„  sin  n9  (41) 

n  2 

and  since  “0,  9  satisfies  the  transcendental  equation 

sin  (N+1)8  «  0  (42) 

Further  since  (41)  is  a  non-trivial  solution,  0  cannot  equal  0  or  tr. 

Thus  rne  roots  of  (42)  are 

. N)  (43) 

n  N+1 


11 


Note  that  they  are  equally  spaced  in  the  opan  interval  (OjiT).  Equation 
(39)  then  yields  the  distinct  eigenvalues 

1-2P  <  X  -  1  +  2p  COS0  <  1+2P  (n-l,...,N) 

n  n 

so  that  is  given  by  (30)  . 

When  the  correlation  of  non-contiguous  pulses  may  not  be  neglected 
we  require  the  eigenvalues  of  the  full  C  matrix.  As  is  easily  verified, 
the  inverse  matrix,  C  is  the  tri-diagonal  matrix 


Thus  the  eigenvalue  problem  to  be  solved  is 


In  the  special  case  N-2,  this  beconres 


which  leads  immediately  to  the  distinct  eigenvalues 


X  -  1  ±  p 


Equation  (30)  then  yields  which  may  be  written  in  the  form 


(1+X)V, 


^  (l^X)  -  <PX)  1  ^ 


OX 


PXV, 


DXV„ 


(1+X)^  -  (DX)^ 


+  cosh 


(1+X)^  -  (DX)^J 


which  is  the  same  result  as  that  of  Schwartz. 


To  solve  (46)  when  N>2,  we  employ  the  tecVinigue  developed  for  the  low 
prf  waveform;  i.e,,  ,we  again  foirmulate  and  solve  ,m  equivalent  homogeneous 
bvp.  Prom  (47)  we  have  the  equation 


-P/ 


pv  -  0  (n-1, - N) 

n+l 


and  the  boundary  conditions 


V  -  pv,  -  0 

0  j. 


v  .  -  pv  •  0 
N+l  N 


Again  there  is  a  solution  of  the  form  (37) ,  provided  that 


_>  1  implies  either  X  2.  or  X  ^  iTp' 

equation  (34)  cannot  be  satisfied.  Thus  we  introduce  the  real  angle  0  by 
means  of 


Since 


1+P 


1-0 

X 


•/ 


(20) 


cos  0 


20 


(55) 


^md  write  the  solution  in  the  form  (40) .  Application  of  the  boundary 
conditions  then  yields  the  pair  of  equations 

[1-0  cosOlKj^  -  (0  sineiK^  •  0  (56) 

(cos  (N+l)6  -  p  cos  N0]K^  +  (sin  (N+l)6  -  p  sin  N01K2  -  0  (57) 

whose  determinant  must  vaunish.  Thus  0  obeys  the  transcendental  equation 
[c.f .  (42) 1 , 


sin(N+l)0  -  2d  sin  NB  +  0^  sin(N-l)6  -  0  (58) 

Since  the  values  9  '0,  ir  do  not  permit  a  non-trivial  solution  to  the  bvp, 
the  roots  of  (58)  again  lie  in  the  open  interval  (0,11)  .  Since  it  has  not 
been  possible  to  solve  (58)  analytically,  the  following  remarks  allow  a 
numerical  solution  to  be  easily  achieved: 

To  show  that  there  are  exactly  H  roots  between  0  and  it  we  first,  write 
(58)  as 


7  2 

((1+0  )cos0-2p]  sin  n6  +  (d-p  )  sin  9]  cos  N6  •  0 

then  introduce  the  function  0(9)  by  means  of 


(59) 


sin  <Ji(6) 


A 


(1-D  )  sin9 


1+P  -  2D  cost 


(60) 


14 


(61) 


cos  (t>(6) 


A 


(l*p^)cos9  -  20 
2 

1  +  p  -  2p  cos@ 


so  that  (59)  becomes 


sin  [n6  +  (f  (0)  ]  -0 

[N.B.  as  p^,  (()(0)  0  and  (62)  (42)] 


(62) 


Since 


d0 


1-P 


1+p  -  2pcos0 


(63) 


we  have 


<t  (IT) 


«(0)  +  / 


1-D 


d0 


1+p  -  2pcosG 


(64) 


Since  (62)  represents  a  modulated  sinusoid  whose  total  phase  increases 
monotonically  from  0  at  6»0  to  (N+l)^  at  ©■IT,  it  has  exactly  N  distinct, 
zero  crossings  in  the  open  interval  (0,Tf)  . 


Further , 


2 

since  ^  (0#^)  and  ^ 

00“ 


1  at  cos6»p,  the  function 


sin  (n6  +((>(6)]  oscillates  more  rapidly  than  does  sin(N+l)9  in  the  domain 
0  <  6  <  cos  and  more  slowly  (but  still  more  rapidly  than  does  sin  n6) 
in  the  domain  cos  <  9  <  tt.  Using  these  observations  concerning  the 
spacing  of  the  roots,  it  is  an  easy  matter  to  accurately  locate  the  roots 
by  means  of  a  Newton-Raphson  method. 

Denoting  the  roots  by  eigenvalues  [c.f.(55)]  are  given  by 


hR 

i+p 


<  X 

n 


RzR. 


l+p  -  2PC086 


l+P 

1-P 


(n«l, . . . ,N) 


(65) 


and  the  detection  probability  by  (30). 

Note  that  expansion  of  (65)  into  a  power  series  in  p  yields  a  first 
order  approximation  which  is  identical  to  (44)  -  a  useful  check. 


111.  RESULTS 


Figures  1-7  present  detection  probability  versus  per  pulse  average 
signal-to-noise  ratio  in  dB  for  false  alarm  probability  .of  10  ^  and  N»2,4, 
6,10,15,20,30  pulses  respectively. 

The  dashed  curve  is  the  Marcum  result  which  was  obtained  by  employing 
(11)  In  (10) ,  performing  the  Laplace  inversion  and  integrating  the  result¬ 
ing  pdf  over  the  threshold;  this  leads  to 


P 


D 


FA 


- (V  +NX) 
e 


OO 


L 

n-1 


n  I 


N-l+n 

E 


k-N 


(66) 


The  closest  dark  curve  is  the  Swerling  case  2  result  (p-0)  which 
was  obtained  by  replacing  by  V.j,/(l+x)  in  (29). 

The  next  (lighter)  curves  are  the  results  for  C“0.40,  0.60,  0.80, 
0.90,  0.95,  0.99  respectively. 

The  last  (dark)  curve  is  the  Swerling  case  1  result  (p-1)  which  was 
obtained  by  performing  the  inverse  Laplace  transform  of  (18)  and  inte¬ 
grating  the  resulting  pdf  over  the  threshold;  this  leads  to 


(67) 


16 


For  N"20,  the  numerical  Instability  exhibited  at  small  values  of  X  when 
0  •  0.99  is  caused  by  the  small  eigenvalues  [c.f.(65)]  which  cluster  together; 
i.e.,  by  indexing  the  eigenvalues  according  to  tize,  it  may  be  seen  that  the 
Individual  terms  in  (30)  alternate  in  sign  and  Increase  in  magnitude  as  the  dif¬ 
ference  becomes  smaller.  This  instability  becomes  more  pronounced  and 

occurs  at  even  smaller  values  of  p  as  N  increases  (it  occurs  at  P  •  0.95  when 
N  “  30).  In  order  to  present  a  clean  figure  when  N  30,  the  curve  for  p  ■  0.99 
has  been  deleted  from  figure  7. 

The  figures  show  that  as  p  increases  from  zero  to  unity,  more  per  pulse  av¬ 
erage  signal- to-r.oise  is  required  to  achieve  the  same  This  increase,  how¬ 

ever,  ii.  smaller  than  one  would  intuitively  expect;  for  example,  increasing  P 
from  zero  to  ope-half  requires  an  additional  Increase  of  less  than  a  dB  in 
signal-to-noise  at  any  0.50  ^  ^  0.99  euid  N  >  1  while,  for  example  at  - 

0.95,  increasing  P  from  zero  to  one  requires  an  increase  of  5.2  dB  for  M  -  2 
or  10.6  dB  for  N  <•  30. 
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F igure  1 . 


'J  -  4  PULSES,  PFfl  "  10mh-6 

NCN-FlUCTUHTINC  TRRGET  RiND 
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Figure  2. 


O  CT) 


N  -  6  PULSES,  PER  -  10mk-6 

NON-FLUCTUflTING  TARGET  Ai\C 
0  -  0,  .4,  .6,  .8,  -9,  .95,  .99,  l.C 


6.0  8.C  10. 0  ;2.0  M.O  16.0  18.0  20. j  22.0  24.0  25.0 

IGiMRL  TU  NOISE  RATIO  IN  DB 


Figure  3. 


N  -  10  PULSES,  PFfl  - 
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PROBABILITY  OF  DCTFCTION 


N  -  15  PULSES,  PER  -  iO^H-b 
NON-FLUCTURTING  TARGET  RND 
RHO  “0,  .4,  .6,  .8,  .3,  .35,  .39,  1.0 
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PART  TWO.  FLUCTUATION  LOSS 


1.  INTRODUCTION 


The  exact  fluctuation  loss  for  a  Gauss-Markov  signal  is  determined  as  a 
function  of  the  number  of  integrated  pulses,  the  correlation  between  conti¬ 
guous  pulses  and  the  specified  detection  and  false  alarm  probabilities.  This 
exact  loss  is  compared  to  Barton's^  approximation  which  calculates  fluctuation 
lota  by  assuming  a  Swerling  case  2  target  and  a  reduced  number  of  "independent” 
pulses.  Two  expressions  for  the  number  (N^  or  N^)  of  equivalent  indepe- -’snt 
pulses  are  derived.  The  domain  of  validity  of  the  approximations  is  estab¬ 
lished. 


See  references  at  end  of  text. 
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11.  BACKGROUND 


Figures  1-7  of  part  I  show  that  at  fixed  and  N  a  greater  per  pulse  aver¬ 
age  aignal-to-noise  ratio  (detectability  factor)  is  required  to  achieve  a  speci¬ 
fied  0.50)  when  the  target  fluctuates  th£u:\  when  it  does  not.  This  re¬ 

quired  increase  ir.  detectability  factor  is  called  the  fluctuation  loss.  For 

fixed  P  and  P  it  depends  on  the  number  of  detected  pulses  and  ♦-.heir  correl- 
D  FA 

ation  and  will  be  denoted  by  L^(M,p). 

Table  I  presents  the  detectability  factors  (N)  for  a  non-fluctuating  tar¬ 
get,  0^ (N)  for  a  Swerling  case  2  (p«0)  target  and  (N)  for  a  Swerling  case  1 
(0-1)  target  at  various  P_'s  of  interest;  also  shown  is  the  normalized  thresh- 
old  voltage  required  to  achieve  the  common  P^^  of  10  .  This  table  has  been 

calculated  from  eqs.  (66),  (29)  with  replaced  by  V^/(l+X).-  and  (67)  of  part  J. 

For  the  P„ ‘ s  of  interest  the  fact  that  for  N>1, 

D  .  . 

L^(N,1)  =  D^(N)  -  Dq(N)  >  D2(N)  -  Dq(N)  -  L^(N,0) 

is  well  knowr.  to  radar  systems  engineers  from  the  previously  published  work  of 
4Marcim\  and  Swerling. 

Whenever  a  surveillance  radar  encounters  the  situation  p«l,  the  target  may 
lie  in  a  deep  fade  for  all  N  pulses  received  on  a  given  scan,  thus  causing  a  de¬ 
tection  to  be  missed.  In  order  to  avoid  this  situation  some  radar  systems  de- 
correlate  the  N  signal  returns  by  transmitting  pulse-to-pulse  frequency  diverse 
waveforms.  This  certainly  decreases  the  fluctuation  loss  from  L^(N,0)  to  L^(N,0) 
(a  considerable  reduction  at  high  and  N)  but  involves  increased  system  com¬ 
plexity  and  cost. 

Fluctuation  loss  depends  not  only  on  P^^  and  P^^  but  also  on  the  fluctuation 
model  and  the  correlation  model.  For  P^^  -  10  a  Gaussian  fluctuation  model 
(Rayleigh  envelope) ,  and  a  first  order  Markov  correlation  model  the  analysis  of 
part  1  yields  the  appropriate  detectability  factors  for  each  N  and  p.  By  si;ib- 
tracting  the  detectability  factor  (N) ,  we  obtain  the  fluctuation  losses, 

(N,  ')  ,  of  table  II. 

It  is  observed  that  the  fluctuation  loss  increases  monotonically  with  P^ 
and,  except  for  the  singular  case  p  -  1,  decreases  monotonically  with  N.  Of 
greater  significance  however  is  the  nature  of  the  monotonic  increase  in  fluc¬ 
tuation  loss  with  D  and  in  particular  its  rate  of  increase  as  p  approaches  1. 
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TABLE  I 


DETECTABLILITY  FACTORS  AT  P 


FA 


10 


N 

?jj-0.50 

Pjj-0.70 

P|j-0.90 

Pjj-0.95 

Pjj-0.99 

V_(dB) 

T 

DqCN) 

DqW 

Do(m 

Dq(N) 

Do(N) 

1 

11.24 

12.09 

13.18 

13.66 

14.50 

11.40 

2 

8.80 

9.61 

10.65 

11.11 

11.91 

12.22 

4 

6.49 

7.26 

8.25 

8.68 

9.44 

13.29 

6 

5.21 

5.94 

6.90 

7.32 

8.06 

14.05 

10 

3.65 

4.35 

5.27 

5.67 

6.38 

15.15 

15 

~  47 

3.14 

4.03 

4.41 

5.10 

16.13 

20 

1.66 

2.31 

3.17 

3.55 

4.21 

16.89 

30 

0.54 

1.18 

2.00 

2.37 

3.01 

18.03 

1 

1 

D2(N) 

D2(N) 

D^CN) 

Dj  (N) 

D2(N) 

1 

12.77 

15.77 

21.14 

31.  S" 

2 

9.52 

11.53 

14.83 

20.47 

i  4 

6.83 

8.28 

10.51 

11.65 

13.97 

1  6 

F  4 

6.65 

8.49 

9.41 

11.22 

'  10 

1 

3.7 

4.80 

6.29 

7.02 

8.40 

,  15 

2.Sv 

3.46 

4.75 

5.36 

6.52 

j  20 

1.71 

2.55 

3.73 

4.29 

5.32 

'  30 

0.58 

1.34 

2.39 

2.88 

3.78 

D^(N) 

Di(N) 

(N) 

D^(N) 

Dj_(N) 

1 

12.77 

15.77 

21.14 

24.29 

31.37 

2 

10.33 

13.32 

18.69 

21.84 

28.91 

4 

8.03 

11.00 

1 

16.36 

19.51 

26.60 

'  6 

6.74 

9.71 

15.07 

18.21 

25,30 

;  10 

5.19 

8.15 

13.50 

16.64 

23.73 

1 

15 

4.00 

6.96 

12.31 

15.45 

22.54 

20 

1 

3.19 

6.14 

11.49 

14.63 

21.72 

1  30 

2.08 

5.03 

10.37 

13. -i 

20.57 

TABLE  II 


L^(N,0)  -  FLUCTUATION  LOSS  (dB) 


N 

p*0.00 

P-0.40 

p-0 . 60 

P-0.80 

p-0. 90 

P-0.95 

p-0. 99 

P-l.O 

2 

0.72 

in 

1.04 

1.28 

mm 

1.48 

1.51 

1.53 

4 

0.34 

WSM 

0.65 

1.00 

1.40 

1.49 

1.54 

6 

0.21 

0.31 

0.47 

0.81 

1.32 

1.47 

1.53 

0.50 

10 

0.12 

0.19 

0.29 

0.58 

0.92 

1.19 

1.43 

1.54 

15 

0.07 

0.11 

0.19 

0.42 

0.73 

1.04 

1.39 

1.53 

20 

0.05 

0.08 

0.14 

0.32 

0.60 

0,93 

1.34 

1.53 

30 

0.04 

0.06 

0.10 

0.22 

0.44 

0.65 

1.26 

1.54 

2 

1.92 

2.13 

2.44 

2.95 

3.31 

3.64 

3.71 

4 

1.02 

1.24 

1.58 

2.26 

2.86 

3.57 

3.74 

6 

0.71 

0.90 

1.21 

1.86 

2.54 

3.53 

3.77 

0.70 

10 

0.45 

0.59 

0.83 

1.40 

2.08 

2.72 

3.44 

3.00 

15 

0.32 

0.42 

0.61 

1.08 

1.71 

2.35 

3.35 

■3.82 

20 

0.24 

0.33 

0.48 

0.88 

1.45 

2.12 

3.23 

3.83 

30 

0.16 

0.22 

0.33 

0.64 

1.12 

1.74 

3. 04 

3.85 

2 

4.18 

4.49 

4.93 

5.79 

6.59 

7.21 

7.88 

8.04 

4 

2.26 

2.59 

3.13 

4.23 

5.36 

6.36 

7.73 

a.ii 

6 

1.59 

1.90 

2.39 

.3.48 

4.66 

5.79 

7.65 

8.17 

0.90 

10 

1.02 

1.27 

1.68 

2.65 

3.78 

5.00 

7.44 

8.23 

IS 

0.72 

0.91 

1.25 

2.11 

3.15 

4.33 

•7.24 

8.28 

20 

0.56 

0.72 

1.01 

1.73 

2.72 

3. 87 

6.98 

8.32 

jH  ' 

0.39 

0.52 

0.73 

1.32 

2, 17 

3.23 

6.58 

8.37 

2 

5.84 

6.33 

7.33 

8.32 

9.18 

10.52 

10.73 

4 

3.36 

3.96 

5.24 

6.60 

7,86 

10,32 

10.83 

0 

2.44 

3.03 

4.28 

5.63 

7.08 

10.22 

10.89 

0.95 

10 

1.64 

2.13 

3.27 

4.62 

6.06 

9.93 

10.97 

15 

0.95 

1.19 

1.60 

2.58 

3.84 

5.26 

9.66 

11.04 

20 

0.74 

0.94 

1.29 

2.16 

3.33 

4.69 

9.32 

11.08 

30 

0.51 

0.67 

0.94 

1.65 

2.66 

3.92 

8.79 

11.14 

2 

8.56 

8.91 

9.47 

10.61 

11.83 

13.04 

16.70 

17.00 

4 

4.53 

4.98 

5.70 

7.22 

1  8.88 

10.56 

17.16 

6 

3.16 

3.60 

4.32 

5.84 

7.56 

9.34 

17.24 

0.99 

10 

2.02 

2.40 

3.02 

4.44 

6.13 

7.93 

15.77 

17.35 

15 

1.42 

1,74 

2.27 

3.52 

5.11 

6.87 

15.34 

17.44 

20 

1.38 

1.86 

3.02 

4.46 

6.16 

14.80 

17.51 

' 

30 

■■i 

0.99 

1.38 

2.31 

3.61 

5,10 

13.94 

17.56 

For  example,’  suppose  all  we  know  concerning  the  target  returns  is  that  they  are 
not  independent;  if  we  make  the  usual  assumption  that  they  are  fully  correlated 
(i.e.,  p»l)  and  integrate  30  pulses  to  achieve  required  P^'s  of  say  C.90,  0.95, 
0.99  respectively,  then  a  system  which  employs  pulse-to-pulse  frequency  divers¬ 
ity  could  recover  the  fluctuation  losses  (L^(30,l)  -  L^(30,0)]  of  respectively, 
8.0,  10.6,  and  17.8  dB.  If  however  p  is  not  in  fact  equal  to  1.0  but  is  instead 
0.9  then  the  recoveradDle  fluctuation  losses  are  only  1.8,  2.2,  2.8  dB  respect¬ 
ively  (at  p-0.8  the  recoverable  losses  are  even  smaller,  viz.,  0.9,  1.1  and 
1.5  dB) .  Thus  the  trade-off  of  system  complexity  versus  recoverable  fluctua¬ 
tion  loss  is  markedly  affected  by  knowledge  of  the  dependence  of  fluctuation 
loss  on  signal  correlation  -  a  dependence  established  in  part  1  of  this  paper. 

In  summary,  the  usually  employed  assimption,  p-1,  should  not  be  facilely 
invoked  simply  because  the  signal  returns  are  known  to  lack  tc.tal  independence. 
Doing  so  results  in  an  extremely  optimistic  estimate  of  the  b'^nefits  which 
can  be  achieved  by  employing  frequency  diversity. 

III.  ANALYSIS 

e  next  consider  two  approximations  which  greatly  simplify  the  determin- 

*7 

ation  of  fluctuation  loss  for  any  correlation  model.  Barton'  introduces  the 
concept  of  the  number  of  equivalent  "independent"  pulses;  this  is  defined  by 

"  Min,  ^N,  It  ^  J  (67) 

in  (67)  is  the  observation  time  and  t^  is  the  correlation  time  (which  Barton 
defines  for  any  correlation  model  as  the  reciprocal  of  the  effective  noise  band¬ 
width  of  the  two  sided  fluctuation  spectrum).  The  definition,  (67)  is  assumed 
to  be  valid  independent  of  or  Barton's  two  conjectures  state; 

Lj(N,P)  L^(N^,0)  (FIRST  CONJECTURE)  (68) 

L,  (N  ,0)  L, (l,p)  (SECOND  CONJECTURE)  (69) 

I  €  N  t 

e 

The  validity  of  these  conjectures  will  now  be  tested  against  the  exact 
results  available  from  our  analysis  of  the  first  order  Markov  correlation  model. 
We  begin  by  calculating  Barton's  N^. 

In  (67)  we  interpret  t^  as  the  correlation  time  of  the  process  at  the  detec- 
^Soc  references  at  end  of  text. 
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tor  cjutput.  Thus  letting  k  and  y  denote  i.i.d.  stationary  Gaussians  which  rep 
resent  the  inphase  and  quadrature  components  of  signal  at  the  detector  input, 
the  detector  output  becomes 


z (c)  -  (t)  +  (t) 

Its  average  (denoted  by  an  overbar)  is 

z(t)  -  x^(t)  +  y^(t)  -  2R  (0) 

X 


(7 


(7 


where  [c.f.(31)] 


Rjj(T)  -  R^(0)  e 


-V  T 


is  the  autocorrelation  of  the  inphase  (or  quadrature)  Markov  process. 
The  covariance  function  at  the  detector  output  is 


R2~(T)  -  [z(t)-z(t)] [z(t+T)-z(t+T))  -  2[P 

X 


which,  since  x(t)  is  a  Gaussian  process,  yields 


(7 


(7 


The  "two  sided  fluctuation  spectrum"  is  interpreted  as  the  Fourier  trans¬ 
form  of  the  covariance  function,  i.e.. 


S„(f)  -  f  4R  ^(0)e 

A  ef  X 


.2vlTl-i2:TfT^^  _  2 

X  .  2 


4v 


C 


4V  +(2TTf) 


This  has  an  effective  noise  bandwidth,  f  ,  given  by 


the  correlation 


Since  the  observation  time  is  (N-l)T  and  [c.f.(31)]  p  -  e 

time  t  ,  is  T/lrt^.  Thus 
c'  '  p 

N  (N,p)  -  Min.  (N,  l+(N-l)ln 

A  more  common  definition  of  the  number  of  independent  integrated  pulses 
is  the  increase  in  the  ratio  of  squared  mean  to  variance  when  N  pulses  are 
summed,  i.e.,  define 


It  is  easily  shown  that 


V 

o  ,  _ 


\  N  t-  2  53  ^N-m)  R^^{in)/R^  (0) 

m“l 


Since  for  the  first  order  Markov  model 


eq.  (80)  yields 


R^(m)  ■ 


Nj.(N,p)  - 


hL  A.I 

1-P^  ^  “  1-0^  ) 


Note  that  Ng  ^  except  atN«lorp«l,  or  P-0  where  equality  obtains. 

Table  III  compares  N^  with  for  the  N  and  p  of  table  II. 

Figures  8-12  are  plots  of  the  exact  fluctuation  losses  of  tcible  II  versus  N^, 
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tha  number  of  equivalent  independent  pulses  as  defined  by  Barton.  The  solid 
curves  represent  L^(N^,0),  i.a.,  the  left  side  of (69)  .  Aside  from  a  neglig¬ 
ible  positive  curvature  which  is  moat  pronounced  at  low  N  euid  high  P  ,  these 

^  o 

do  (as  Barton  conjectures)  exhibit  a  linear  decrease  with  N  .  The  dashed  curves 

represent  the  right  side  of  (69),  i.e.,  —  L  (l,p) .  We  see  that  the  solid  and 

N  r 

e 

dashed  curves  are  within  1/2  dB  of  each  other;  further  we  note  that 


<  ^ 
N 

e 

L^(1,P) 

(Pp  -  0.50) 

(83) 

L, (N^,0) 

I  Q 

>  2^ 
N 

h^(l,0) 

(Pj^  -  0.70,  0.90,  0.95,  0.99) 

(84) 

e 


In  Figures  10-12.  the  upper  solid  curve  is  the  continuation  of  L^(N^,0) 

into  the  region  10  <  N  ^30  (divide  the  ordinate  by  10  and  multiply  the 
~  *  1 

abscissa  by  10).  The  dashed  curve,  —  L^(1,0),  is  its  own  continuation. 

The  symbols  •,  x,0.  A,  O  ,  V,  identify  values  of  L^(N,p)  for  p  ■  0.40,  0.60, 
0.80,  0.90,  0.95,  0.99  respectively. 

we  see  that  for  ■  O.SO  and  0.70,  Barton's  first  conjecture  underesti¬ 
mates  tlie  exact  fluctuation  loss  by  a  fraction  of  a  dB.  For  "  0.90,  0.95, 
0.99  the  first  conjecture  overestimates  the  exact  fluctuation  loss  at  small 
(by  more  than  a  dB  at  •  0.95  and  by  several  dB  at  P^  -  0.99)  and  xmder- 
estimates  the  exact  fluctuation  loss  by  less  them  a  dB  at  large  N^.  Since 
the  fluctuation  loss  monotonically  decreases  with  N^,  Barton's  first  conjec¬ 
ture  gives  a  result  which  is  conservative  when  the  losses  are  large,  and  op¬ 
timistic  when  the  losses  are  small. 

If  one  were  to  use  N  instead  of  N  to  represent  the  number  of  indepen- 
1  e 

dent  pulses,  the  symbols  representing  L^(N,c)  would  shift  to  the  left  so  that 
the  approximations  would  become  even  more  conservative  at  *  0.90,  0.95, 
0.99  when  the  losses  are  large.  Thus  Barton's  is  the  preferred  approxi¬ 
mation  for  0.90  (large  losses),  while  the  approximation  using  ir.  to 

be  preferred  for  ^  0.70  (small  losses). 
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NUMBER  OK  INDEPENDENT  PULSES 


Figure  9. 
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SUMBFR  OF  IMDEPESDEST  PULSES 


Figure  10. 
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